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FOREWORD 


This document is tlie first in a five volume report which describes a 
comprehensive digital computer simulation of the dynamics of heavy lift 
airships and generically similar vehicles. 

The work. v;as performed by Systems Technology, Inc., Hawthorne, Cali- 
fornia for the Aeronautical Systems Branch in the Helicopter and Powered 
Lift Division of the National Aeronautics and Space Administration, Ames 
Research Center, Moffett Field, California. The simulation development 
was carried on between September 1979 and January 1982 and is currently 
installed on the Ames Research Center CDC 7600 computer. Je contract 
technical monitors for NASA were Dr. Mark Ardema, Mr. Ai.sn Faye, and 
Mr. Peter Talbot. STI's Program Manager was Mr. Irving Ashkenas. 

The authors wish to acknowledge the technical contributions of 
Mr. Robert Heffley, Mr. Thomas Myers, and Mr. Samuel Craig and the fur- 
ther contributions of Mr. Aliya Hall, Ms. Natalie Hokama and Ms. Leslie 
Hokama in simulation software development. Special thanks are due to 
Ms. Kay Wade, Ms. Linda Huffman, Mr. Charles Reaber, and STI's produc- 
tion department for the preparation of tlia five volumes of this report. 
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BACKBIOaMD AMD ODJICTITIS 

Tht h««vy lift airship (HLA) la a vahlcla (Fig. 1) In vhlch tha 
tachnologlas of rotorcraft and Ilghtar than air (LTA) art coablnad to 
provlda a aaana of lifting, transporting, and placing axtraaaly Basslva 
payloada ■— up to parhapa 100 tons* Tha lifting powar of tha buoyant 
gaa of fasts tha wslght of tha structura and propulsion syatams; tha 
lifting powar of tha rotors supports tha payload and provides vahlcla 
control. Tha concept has racalvad consldarabla attention In racant 
years with technology ravlavs and economic feasibility studies, 
Refs. 1-8. 

Tha operation of these vehicles In their Intended roles, and the 
associated engineering problesM, have been studied to a much more 
Halted extent, e.g., Refs. 9 and 10. Like the classic airships, tha 
response to gusty environments la an Important Issue which has received 
Halted analytical treatment (Refs. 11-13). Unlike Its pradecesaors, 
Che HLA will typically be used In missions requiring prc.ls* control 
through flight regimes encompassing large and rapid changes of speed. 
Incidence, and load. There la a need, therefore, for an In-depth under- 
atandlng of the flight mechanics, handling qualities and flight control 
requirements of the HLA. 

The primary objective of Che technical effort described In this 
report la Che development of a generic, yet comprehensive mathematical 
model and computer simulation of the HLA flight dynamics over Its antlre 
flight envelope. ImpHclC In this simulation development are the data 
reviews and analyssf which suppor' tha equations of motion and the 
calculation of forces and momenta acting on the vehicle. The slmula- 
tlcn, HYBRDS, it addressed to the broad requirement noted above and la 
Intended for use as a aynthasla and analysis tool for tha evaluation of 
competing design concepts. 
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Figure 1* Tjrplcal Heavy Lift Airship With 
Slung Payload 


StMMAV.T SUfOLATIOH DESCRIPTION 

The HYBROS (Hybrid Buoyant Rotorcraft Dynamic Simulation) Is a set 
of computer programs that perform non- real- time, nonlinear simulations 
of the flight dynamics of HLA-llke vehicles — hybrid multi-rotor air- 
craft* The three main programs available with HYBRDS are: 

• HLASIM — models the powered vehicle in flight. 

• HLAPAY — models the powered vehicle In flight 
cari.ylng a slung payload. 

• HLAMOR — models the unpowered vehicle moored to a 
rigid support at a single point. 

Vehicle configurations that can be simulated Include free balloons, 
classical airships, hybrid airships (HLAs), and non-buoyant multi-rotor 
helicopters. The modeling Includes the so-called apparent mass effects 
of the buoyant hull, a number of Interference effects arising out of the 
proximity of hull, rotors, tail surfaces, and ground plane, and quasl- 
statlc rotor coring and flapping. The formulation of the equations of 
motion allows calculation of loads between the main Interfacing parts of 
the vehicle structure. 
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The basic configuration la the buoyant quad- rotor concept — four 
llft-propulalon unlta (LPUa) dlapoaed about a central buoyant hull which 
has controllable tall aurfaces* Each LFD has a lifting rotor and a 
thrusting rotor (propeller), both with collective controls; the rotor 
has lateral and longitudinal cyclic controls as well* Other configura- 
tions are achieved by selectively "turning off” various portions of the 
model by altering the Input data* 

The HYBRDS programs each perform three major functions: 

1) Calculation of operating point or trim condition, 
e*g*, for performance estimates. 

2) Calculation of stability and response derivatives 
about selected operating points for use In linear- 
ized vehicle dynamic response models. 

3) Forward Integration la time of the vehicle state 
variables In response to control commands and 
aerodynamic disturbances, starting at selected 
trim conditions. 

In addition to the independent state variables, the output data 1 '-'.lude 
a large number of dependent variables such as aerodynamic (and other) 
load contributions, rotor thrust levels, and control surface deflec- 
tions * 

The simulation was developed on a Cyber 176 using the Interactive 
NOS operating system* It was written In ANSI FORTRAN IV and contains 
some 65,500 statements In 343 subroutines, most of which are shared 
among the three main programs. The simulation requires the IMSL library 
for numerical processing support. It is operational on CDC 7600 com- 
puter under the SCOPE 2*1 operating system at the Ames Research Center 
of the NASA* The source code Is heavily commented to facilitate addi- 
tions and modifications by the using organization to suit configuration- 
peculiar requirements* 

DOCUMENTATIOH 

The simulation program documentation Is contained In the final 
report, "Flight Dynamics Analysis and Simulation of Heavy Lift 
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Airships*" The report Is in five volumes, of which this volume Is th. 
first* The volumes carry subtitles Indicating the Intended readership: 

Volume 1 Executive Summary 

Volume II Technical Manual 

Volume III User's Manual 

Volume IV User's Manual, Appendices 

Volume V Programmer's Manual 

The following paragraphs briefly describe the contents of each* 

This volume, the Executive Summary , provides an Introduction to the 
remaining volumes* It contains a description of the simulation In 
sufficient detail to allow the reader to determine the potential appli- 
cability of HYBROS to his own use* It describes key features of the 
technical development and outlines potential additions i-nd modifications 
to the program foreseen at the time of writing* 

The Technical Manual describes the mathematical models contained in 
the simulation in considerable detail with supporting evidence for the 
model forms chosen* It describes the trimming and linearization algo- 
rithms used in the simulation* ^pendlces to the manual outline the 
coefficient estimation procedure for the input data and provide example 
simulation results* 

The User's Manual provides the basic information necessary to run 
the programs* This Includes descriptions of che various data files 
necessary for the program, the various outputs from the program and the 
options available to the user when executing the program* Additional 
data file information Is contained In the three appendices to the 
manual* These appendices list all Input variables and their permissible 
values, an example listing of these variables, and all output variables 
available to the user* 

The Proara mer's Manual is Intended for the maintenance programmer 
who will support the program* It contains explanations of the logic 
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rtabodi«d in th« various prograa aodulas* a dictionary of program vari-' 
ablrs, a aubrouclna liating, aubrouclna/coar •n-block/crota'-'rafarartc*'. 
listing, and a calling/callitd subroutlna croaa-rafaranca listing, fha 
■anual doss not rapaat data aval labia in tha Uaar*a Manual . 

Wind tunnal taats of rha quad-rotor daaign concapt wara originally 
planned to cotneida with tha tlaa span of tha siaulation davalopsMnt. 
Tha purpose was to provide data currently lacking in tha litaratura, 
particularly for hull-rotor intarfaranca affects. This effort was tar- 
alnatad before the wind tunnal nodal was built. Tha daaign considera- 
tions and modal plans are contained in two technical sMaoranda, Refs. 14 
and IS. 

In addition, tha formulation of tha nathanatical aodals, tha soft- 
ware tasting, and simulation axarcisa resulted in t number of papers, 
Rafs. Ib-19. The subject of most of these papers was a hypothetical 
four rotor, four thrusting propeller riLA configuration having a hull of 
low fineness ratio and small tail surface area. This particular con- 
figuration was defined early in the program to serve as a test case 
throughout simulation development. 
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OVERVIEW 

The heavy lift airship can be described as resembling a clasialc 
airship design where the propulsion scheme has been extended and elab- 
orated* The several engines of the airship are replaced with consider- 
ably larger lift-propulsion units (LPUs) mounted on relatively long 
outrlgger-llke arms to provide clearance for the main rotor of each LPU* 
The LPOs also carry a propeller to provide propulsive force and/or addi- 
tional control force, particularly when the main rotor Is very lightly 
loaded* The payload. If any. Is suspended below the HLA with a multi- 
cable sling system* 

The most complex parts of the simulation are associated with the 
aerodynamic modeling of the vehicle* In calculating the aerodynamic 
loads acting on various parts of the vehicle structure, the emphasis Is 
placed on determining the dominant effects using aerodynamic models 
based on uniformly valid first approximations* The details of the 
models are described fully in Volume II of this report; the Intent here 
Is to provide a broad catalog and explanation of the major features* 

The simulation also Includes models of the suspension C6> the 

landing gear, and the flight control system* In all ses t’ *pe are 
quite simple in keeping with the overall philosophy of exprsing the 
dominant properties of the generic HLA concept* The source code for the 
flight control system Is confined to a few subroutines to make changing 
the system relatively easy; the flight control system Is not "universal" 
because of the broad range of possible HLA control configurations* 

MDLTIBOOT EQOATIOHS OF NOTION 

In the simulation equations of motion, the HLA Is treated as an 
assembly of Interconnected rigid bodies — the central hull and support 
structure assembly, and the several attached LPUs* The mathematical 
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represencaclon of similar structures has received considerable attention 
In recent years, particularly in connection with spacecraft dynamics, 
e>g>, Refs. 20-22. The objective in most of these investigations Is to 
arrive at a computationally economical set of equations for a specific 
configuration within a general class of related arrangements. 

In the case of the heavy lift airship, the principal requirements 
are to provide for variability In the configuration and to give some 
measure of the loads between components. For these reasons each body Is 
treated as a separate, albeit constrained, entity with Its own axis 
system and equations of motion. There are additional constraint equa- 
tions which express the relationships between the accelerations of pairs 
of bodies resulting from their interconnections. 

To the user, this treatment of the rigid body modeling of the hull 
and attached LPUs is reflected in the input data requirements which call 
for the geometrical and inertial data of each body — the hull and each 
of the attached LPUs — rather than data for the structure as a whole. 
This facilitates changes In, for example, the location and orientation 
of the LPUs In predesign studies. One need n> recompute the Inertial 
properties for the structure as a whole. Th.s approach requires r 
defined net center of gravity for the entire HLA; the constraint for( «i 
and moments between the hull and each of the LPUs cause the structure to 
move as a rigid body. 

When solved simultaneously in the program, the equations of motion 
and the equations of constraint yield solutions for the vehicle accel- 
erations in the independent degrees of freedom as well as solutions for 
the forces and moments acting between pairs of bodies which provide the 
constraint. When the Interfaces between the LPU and hull-plus-support 
(outrigger) structure are appropriately chosen, the result is a readout 
of these internal loads as well as motions of the vehicle. 

AERODYNAMIC CALCULATIONS 

The aerodynamic modeling of the HLA is similarly organized. Forces 
and moments are calculated for the hull (tall off), tall assembly, each 
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LPU, and Che payload. If present* Figure 2 shows the conputaClonal flow 
for these calculatloii^r 'la various steps of which are as follows: 

1} The vehicle states are taken from the current time step; 
these Include translational and rotary motions of each 
element, local alrmass motions at each element (steady and 
turbulent wind Inputs), and several control Inputs* 

2, 3) These are resolved Into relative air-to-element velocity 
components at each rotor and propeller hub, each LPU aero- 
dynamic center, the hull cenfr of volume, the tall refer- 
ence center (the aerodynam. . center of effective tall- 
plus-fuselage ensemble) , and the slung payload. If 
present* 

4) The relative velocities are adjusted by various factors or 
Increments to account for Interference due to: 

a) hull on rotors and propellers 

b) rotors and propellers on hull and tall 

c) rotors on propellers and LPU nacelles 

d) ground plane proximity to hull, tall, rotors and pro- 
pellers 

Thus the net relative or apparent velocities are the 
vector sum: 

yrelatlve alrmass 

5) The aerodynamic forces and moments are computed as func- 
tions of the relative velocity at each element* Inter- 
ference effects that arise from changes In the nature of 
the local flow (e*g*, rotor- Induced turbulence In hull 
local flow) are accounted for In the equations of the 
respective element* Buoyant forces are computed at the 
hull center of volume from the normal atmospheric pressure 
gradient, and the horizontal pressure gradients arising 
from changing wind velocities or convergent wind fields. 

6) Net hull forcen at the center of volume (c*v*) are summed 
and transferred to the hull center of gravity (c*g*) along 
with the tall-on-hull forces* At this stage the major 
pitch and yaw stability effects of the hull and fin assem- 
bly are apparent* 

7) The rotor and propeller forces and moments are transferred 
to the LPU c*g* for use In the muxclbody equations of 
motion* 

8) The various loads and control deflections are Inputs to 
other parts of the program* 
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Figure 2. Computational Flow of Aerodynamic Subroutines 
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HULL AND TAIL AEKODYNAHICS 


Quasi-steady and unsteady aerodynamic loads are accounted for in the 
simulation for flight conditions encompassing large ranges in speed » 
incidence angle, and turbulence level. The modeling uses simple func- 
tional dependencies on body axis linear and angular velocities and 
accelerations relative to the alrmass. The dependencies were selected 
to allow approximation of typical empirical data rather than being esti- 
mated from vehicle geometry. In the case of the tall surfaces, pre- 
stall, transition and post-stall flow regimes are included. 

The quasi-steady aerodynamic models were validated against extensive 
wind tunnel data for certain of the classic airships, Refs. 23-27. Fig- 
ure 3 shows an example of this comparison with data for the airship 
Akron for incidence angles out to twenty degrees. The closest match is 
for the drag data which s within 5 percent. Lift is matched within 
10 percent for the hull-plus-tail configuration. The pitching moment 
data require adjustments (X and n) to capture the low incidence insta- 
bility and high incidence stability so typical of airships. Additional 
correlations with high incidence and oscillatory (damping) data indicate 
that the hull and tail aerodynamic forces are generally valid to within 
about 25 percent over the entire HLA operational envelope. Since nearly 
all of the aerodynamic problems requiring control are simulated in the 
model, further complexity was deemed unwarranted. 

T!ie unsteady aerodynamic forces are those arising out of vehicle 
acceleration with respect to the airmass and accelerated airflow due to 
gusts and irbulence. The effects modeled in the simulation are the 
classic accelerated motion forces (Refs. 28 and 29) which were supplied 
in the analysis of the classic airships (Refs. 30 and 31). Accelerated 
motion of the airmass also leads to horizontal buoyancy forces because 
of the pressure gradient associated with the acceleration. These, too, 
have been accounted for. However, the unsteady lift arising out of 
changes in circulation flow represents a consider bly smaller contribu- 
tion and has been neglected. 
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Figure 3. Simulation Quaal-Steady Aerodynaralce Compared 
vlth "Akron" Wind Tunnel Data 
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The unsteady aerodynamic forces can be a doalnanc contribution Co 
the overall aerodynamic loads on the airship hull* Figure A shows an 
example simulation time history where Che gust is applied only at Che 
tail assembly. The HLA is not actively controlled. The comparison 
between the quasi-steady force contribution (Curve C2) and the total 
force (Curve Cl) shotfs the dominance of the unsteady contribution (Curve 
C3). Note that the gust cannot include step discontinuities in velocity 
as this implies fnflnit.^ acceleration. 

LIFT-PROFULSION UNIT ABROOmNlCS 

The individual lift-propulsion units may be likened to a helicopter 
where the tall rotor has been replaced by a thrusting propeller oriented 
in some arbitrary direction. The aerodynamic forces of Importance are 
those originating in the main rotor and propeller and, to a much lesser 
extent, the LPU nacelle (helicopter fuselage). The control inputs are 
the rotor collective, lateral cyclic, longitudinal cyclic and rotational 
rate; and the propeller blade pitch (collective) and rotational rate. 

The main rotor aerodynamics follow the modeling of Bram 11 in 
Ref. 32. Quasi-steady coning and flapping of the rotor blades is 
assumed with the flapping hinge located at the rotor hub. The inflow 
velocity is uniform across the rotor disk and is assumed quasi-steady. 
An iterative procedure Is used to solve for the Inflow velocity and 
thrust coefficient; torque and flapping coefficients follow directly. 

An unusual added feature Is the replacement of the Bramwell model by 
flat plate drag when the Inflow approaches zero. This approximates the 
vortex ring state which Is frequently encountered on an HLA, particu- 
larly when the vehicle Is not carrying a payload. The modeling is 

not accurate In this region and the program flags this with a message 
telling the user when this flow condition occurs. 

The propeller aerodynamics follow those of the rotor, deleting the 
cyclic controls. The nacelle Is modeled as a bluff body having drag 
along each reference axis dependent upon the relative velocity along 
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Figure 4. Time History of Aerodynamic Terns IXiring a 
One-Mlnus-Coslne Sldeirards Tail Gust 
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that axis. In all cases (rotor, propeller, nacelle), the forces calcu- 
lated depend upon the local velocity relative to the alnaass — at the 
rotor hub, propeller hub, and nacelle aerodynamic reference center. 

PAYLOAD AERODYNAMICS 


The payload aerodynamics are those of a rectangular container, 
neglecting unsteady (Strouhal) flow forces. The aerodynamics are some- 
what more complex than those of the LPU fuselage in that moment terms 
dependent on linear and angular velocites are included. 

XNTEBFEREHCB effects 

The proximity of the rotors and propellers to the hull the ground 
plane, and the tall surfaces leads to mutual interference effects of two 
kinds: those due to velocity changes and those due to the change in 

flow angularity. The former are modeled as velocity decrements; the 
latter by changes in certain of the aerodynamic coefficients. Since the 
modeling here as elsewhere concentrates on first-order effects, it was 
possible to avoid iterative calculations. The functional form of the 
models was formulaL^ 2 d from the available literature on ground effects, 
e.g., Refs. 27 and 33, plus recent flow simulation work. Ref. 3A, and 
ultimately by recourse to first principles (e.g., potential flow solu- 
tions, linear single dependency models, etc.). 

AERODYNAMIC DISTURBAMCES 


The aerodynamic environment is made up of three separate components: 

1) Discrete, one-mlnus-coslne gust inputs acting on 
isolated vehicle elements (e.g., tail only as in 
Fig. 4). These can be applied at several points 
at arbitrary time intervals to represent gust 
"waves." They translate and rotate with the vehi- 
cle and are Intended for test purposes only. 

2} Steady wind which is constant with respect to the 
non-rotating inertial reference frame. The wind 
contributes to the unsteady aerodynamic forces 
acting on the hull and tail when the vehicle is 
rotating with respect to Inertial space. 
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3) Spatially distributed random turbulence which Is 
approximated by a four (or five when payload Is 
present) point source model. 

The aerodynamic force calculations require local relative velocities 
at widely separated points on the vehicle. These locations are far 
enough apart that extrapolation of point gust values and gradients from 
a single point gives Inaccurate results (Ref. 35). The degree of gust 
component correlation between any pair of locations drops off as the 
space between them increases, so that at 100 ft apart — a dimension 
typical of the distance between rotors in a quad-rotor HLA — the corre- 
lation nearly vanishes. 

In the simulation, the hull Is treated as a spatial averager of 
gusts that originate at four points around a meridional plane. They are 
selected by the user to be close to the rotors, the tall surfaces, and 
the reference length of the hull. Appropriately weighted averages of 
these four sources are taken as the gust component *'seen" by the LPUs, 
by the hull, and by the tail assembly. This same model also provides 
gradients along and across the hull which gives rise to rotary gust and 
alrmass acceleration terms. 

The model requires that a set of twelve inertial velocity component 
time histories be computed ahead of time, three mutually perpendicular 
components at each of the four source points. For computational conven- 
ience, the Input points translate but do not rotate with the vehicle. 
This allovjs use of actual gust records, or artificially generated time 
histories containing both correlated and uncorrelated components. The 
user might want correlated components to represent local airflow In 
close proximity to Irregular terrain (e.g., a ridge, near a large 
hangar, etc.). 

The payload. If present, requires an additional six components, 
three translational and three rotational terms, the latter being re- 
quired because of the single-point nature of the payload gust modeling. 
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PATLOAD SLIHG DTMAHICS 


In the HLAPA7 simulation the payload Is suspended below the hull on 
as many as four suspension cables* each modeled as a spring-damper sys- 
tem* The geometry la arbitrary irlthln constraints Imposed by the trim- 
ming routines* Generally this means that sling geometries that do not 
restrain payload angular motion are not allowed and must be approximated 
u> geometries that do provide these restraints* see Volume III* 

This modeling approach was chosen over an alternate rigid link con- 
straint modeling for reasons of simplicity* The rigid link approach 
requires a rework of the multibody equations of motion and constraint 
equations* The constraint equations change for each variation In sling 
topology* e*g.* two versus four cables* blfllar versus pal .red Inverted 
vee* etc** The flexible link approach adopted requires only a single 
computer code configuration to represent arbitrary geometries Involving 
up to four cables* 

LANDIHG GEAR DYNAMICS 

Each of the as many as four landing gears Is modeled as a vertically 
oriented spring-damper combination with two gradients* A Coulomb fric- 
tion model Is used for the forces acting at the ground contact patch* 
The stronger of the two gradients accounts for structural deflection 
when the weaker spring has bottomed out* 

MOORING SYSTEM 

In the HLAMOR simulation, a single mooring attach point Is modeled 
as an additional constraint on the hull motion; the set of constraint 
equations Is augmented* The HLA Is modeled with zero rotor and pro- 
peller speeds; the aerodynamics are replaced by simpler models of the 
at-rest forces acting on the rotor and propeller* 

FLIGHT CONTROL SYSTEM 

The flight control system operates the rotor* propeller, and tall 
surface deflection controls through the software equivalent of a mixer 
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box* The exact nature of the mixing Is alterable by the user by making 
changes within the source code of a single subroutine* 

The mixing "organizes" the many individual vehicle controls into six 
approximately orthogonal linked controls, one for each of the six 
degrees of vehicle motion freedom* Thus changes in the vehicle control 
configuration being analyzed by the user require changes in simulation 
code* 

The mixing is effective not only for operation of the flight control 
system but also for establishing vehicle trim conditions* The trim rou- 
tine operates through this mixing as does the determination of the 
linked control response derivatives* 

The flight control laws themselves are of the classic proportional. 
Integral, derivative (PID) type where the integrator output defines the 
trim deflection of the linked '''^ntrol point* Figure 5 Illustrates the 
operation of this system in regulating against discrete downward gust 
at the tall* The sensed quantities include angular rates, body atti- 
tudes, accelerations and speeds (either inertial or airmass referenced)* 
The system Includes limits on the integrator output, on the linked con- 
trol deflection, and on each of the Individual control points of the 
vehicle* Thus the simulation is capable of inveotlgatlng situations 
where limited control power is a factor* 

The system accommodates user command Inputs at several points* The 
individual control surfaces and the linked control points can have 
finite-duration pulses applied for test purposes* The flight control 
loops themselves can be commanded by user-specified time sequences cf 
commands for execution of mantjuvers, for example* 

The system design is not Intended to be universal; quite the con- 
trary* The generic nature of the EDLA configuration and the range of 
missions for which it is designed both imply control objectives that 
vary over a broad range* As presently defined, the system is oriented 
toward the hover condition with user-specified limits, gains, and sensor 
selection (from a limited set); the software mixer box is Intended for 
four rotor, four propeller configurations* Source code changes in a 
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Figure 5* Conparlson of Open- and Closed-Loop 
Responses to Tall Down Gust 
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SOFTIUIE DBSCRIPTUMi 


SOUBCK a»B mVELOnOBR 

The source code was developed in a top-down fashion wherein the data 
structure. Input/output characteristics, basic program algorithm, etc*, 
were defined early to establish a framework for the development that 
followed* The major modules, each consisting of an Interrelated set of 
subroutines, were developed in sequence, then incorporated and tested 
with the main program(s)* 

The individual modules were developed by a four-step sequence of 
model definition, code development, software test, and comment prepara- 
tion* Model definition required establishing the pertinent equations 
and data requirements* This was often a lengthy process Involving care- 
ful consideration of the relative merits of competing models* Some- 
times, e*g*, hull/tsil aerodynamics and hull/rotor interference effects, 
it required extensive literature review and analysis to establish the 
form of the model and the relevant data structures and ln:erfaces with 
the remainder of the program* 

Reducing the model to FORTRAN source code was straightforward, hut 
the testing that followed often was not* "Stiff" equations, e*g*, land- 
ing gear spring-dampers, resulted in the expected problems with trimming 
the vehicle and Integrating the states forward in time unless the par- 
ameters were chosen with discretion* During and following the testing 
of the individual modules the code was heavily commented by the engineer 
and programmer Involved* These comments complement the program documen- 
tation given in these volumes* 

The development of HLASIM followed the above procedure, the first 
major module being the multibody equations of motion* This was followed 
by the trimming module and the LFU aerodynamics* The modules for the 
hull and tall aerodynamics, the stability and response derivative 
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determination, the flight control systen, and the test gust inputs 
followed In that order* 

Adding the slung payload began with an entire separate program In 
which the dynamics of the slung payload were modelf»d In Isolation — 
payload equations of motion, aerodynamic and cable forces, and trimming 
routines. Following verification. It was "merged" with the HLASIM pro- 
gram to create HLAPAY. This Involved revisions to the trim, stability 
and response derivative and Integration routines. Interference effects 
and movable tall surfaces were then added to both programs. 

HLAMOR was developed by assembling pieces of HLASIM, adding the 
mooring point constraint, and zero speed rotor and propeller aerodyna- 
mics to the system. Again the trim and stability and response deriva- 
tive routines were specialized for this simulation program. 

HYBRDS development was completed by recasting HLASIM and HLAPAY Into 
an overlay structure and adding two small post-processor routines. 
HLAMOR Is small enough to fit on the CDC 7600 without overlays. The 
programs were Installed and checked out on the host computer by compar- 
ing test runs with previously computed check cases. 

MAIN PROGRAM ALGORITHM 

Each of the three programs has the same basic execution algorithm: 

• Read vehicle data and Initialize 

• Read trim specifications 

• Calculate trim state and print 

• Calculate derivatives and print (If requested)—-* 

• Read time history specifications and Initialize 

• Integrate forward one time step — ' 

• Print output at print Interval — — — » 

• Stop 
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The program first reads In Che geometric • Inertial and aerodynamic 
data describing the HLA configuration and Initializes the various data 
arrays. Following this, additional data specifying the trim condition, 
which Includes atmospheric parameters, la read In and a trim state Is 
computed. The printout that follows Includes trim routine performance 
measures to Inform Che user of Sny difficulties encountered In trimming, 
e.g., limiting of the control surfaces. The program now loops back to 
calculate another trim condition. 

Figure 6 shows an example of a trim state sequence, In this case 
Illustrating lateral cyclic deflection requirements as inctloc of 
crosswind '■'eloclty for an HLA with and without slung pa/ At each 

trim condition. If Che user has so requested, the program calculates the 
stability and response derivative matrices associated with the trim con- 
dition and prints them out before looping back to calculate another trim 
state. For each set of derivatives, the program also calculates the 
eigenvalues and eigenvectors of the system. Figure 7 shows an example 
of how certain of the eigenvalues vary with flight condition In a typi- 
cal HLA. 

The structure of Che program to this point supports analyses of 
vehicle performance (Fig. 6) and flight dynamics (Fig. 7). The trim 
specifications can cover the Intended operating envelope of Che HLA; the 
results can point out potential performance or control authority prob- 
lems. The linearized dynamics support studies of flying qualities and 
control system synthesis. Following calculation of the last trim or 
derivative array, the program stops unless the user has requested a 
simulation run* 

If a time history has been requested, the program reads In the 
specifications that Include the flight control system parameters, the 
commands, and Che gust disturbances. The run begins at the last trim 
condition calculated. The program Integrates Che staf.^ vector rate 
of change forward one time step, writes all program output (approxi- 
mately 1000 variables In the largest program, HLAPAY) to a plot file, 
then loops back to calculate another time step. This process continues 
until the print Interval Is encountered at which point user-selected 
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Figure 6« Trimming for Hover In a Steady Crosswind 
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Figure 7. Motion of HLA Complex Roots with Airspeed 

output variables are printed out. The program then enters the integra- 
tion loops once again. The sequence ultimately stops when the maximum 
simulation time specified by the user is encountered. 

TRIM ALGORITHMS 

The trim algorithm determines the elements of a control vector that 
set pertinent elements of the state vector rate of change to *ero. I' 
Is based on a generalized secant method described in Ref. 36. To accom- 
modate wide variations In the HLA mathematical model, it starts with an 
estimate of the needed control vector which is necessarily quite crude. 

Each of the three programs has a different implementation. In 
HLASIM, the control vector '•orresponds to the six linked controls of 
the vehicle. In HLAPAY, the payload is first trimmed in isolation, 
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Che control vector being the three angular components and three linear 
components of payload position. This results in fixed cable forces on 
the hull, which are Chen trimmed with the same six linked controls as 
before. In HLAMOR, the control vector consists of the three angular 
components of hull position. 

The trim condition printout is preceded by measures of Che trim 
algorithm's performance. These data are intended to provide clues to 
the user in the event that the trimmer has difficulty in establishing 
the trim condition. 

STABILITY AND RBSPCNSB DERIVATIVE ALG(»ITHN 

The derivatives are determined by forward and backward differences 
in the state vector rate of change and in a selected set of dependent 
variables, e.g., the constraint forces between the LPUs and the hull. 
The printout of these arrays is followed by an identification of those 
array elements where the forward and backward differences are marked. 
This is intended to flag instances of pronounced nonlinearity around the 
operating point. An IMSL routine is used to calculate the eigenvalues 
and eigenvectors. 

INTEGRATION ALGORITHM 

The integration routine, DVERK, used in the program is based on the 
Runge-Kutta-Verner method. It is an IMSL routine documented in Ref. 37. 
This routine varies the size of the time step within user-specified 
limits to satisfy its error criterion. Operation of the routine there- 
fore requires Interpolation between successive data points defining the 
time history of the commands and disturbances. Even though these may be 
random disturbances, they must be predetermined from the standpoint of 
the integration routine. Gust inputs and control sequences are required 
to be data strings established prior to the beginning of a run, i.e., 
"prerecorded" data. 
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FOST-PtOCKSSOR BOOTUISS 


Two post-processors are included in HYBRDS* One foraats dat<* for 
plotting by utility routines on the Aims 7600 aachine* The other 
aodlfica the gust input response derivative aatrix so chat the effective 
input points are the four or five gust sources rather than six points on 
the hulI-LPU asseably and a seventh on the payload* 

IKPOT 

The input file structure is such that a coaaon set of flics accosiao- 
dates all three prograaa within HYBRDS* Which of the three prograas is 
run is specified In an operating system procedure file* Instructions 
for the programs are contained within a file named INPUT* This file 
specifies the number of trla conditions to be calculated, whether or not 
a time history is to be calculated, and the system of units — English 
or metric* 

The vehicle data are contained in GHDTA (geometry and inertial 
data), ARODTA (aerodynamics of hull-LPU assembly), IFCDTA (Interference 
effects data), and PLNDTA (control surface limits and rotor/propeller 
speeds)* The trim states are calculated using the specifications con- 
tained in TRHDTA* This file also flags whether or not derivatives are 
to be calculated for the Individual states* 

Time histories are calculated per the flight control system and com- 
mand specifications in HISDTA, and the gust input strings for the four 
gust sources contained in four files, RG1-RG4* The printout variables 
are specified by code numbers contained in OUTLST* 

The foregoing is applicable to running -:he HLASIK program* To run 
HLAPAY, additional files are needed* PAYDTA contains the payload geome- 
tric, inertial, and suspension cable dat< RG5 and RG6 contain the pay- 
load translational and rotational gust components and PYOUTL identifies 
the additional, payload-peculiar variables desired in the output list- 
ing* To run HLAMOR, the data file MORDTA is needed, which contains the 
mooring trim specifications and initial offsets from the trim condition 
which are used to excite the mooring dynamics* 
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If th* u««r h«« to rc^u««t«d (In INPUT), all of th« output vurlablo* 
«r« wrltton to binary forsat flla, PLOT, aftar aach trla calculation an^ 
aach intagratien tiaa atap> 

Tha prlncad output atarta with a run daacription — tha coaaanta tha 
uaar haa inaarti«d in INPUT* Following thta ia a labalad liating of all 
tha Input variablaa (axcapt tha guat input atrlnga) if the uaar haa ao 
raquaatad* 

Naxt ara tha raaulta of tha trla calculationa — tha variablaa 
apaclflad In OUTLST and PTOUTL for aach tria condition, togachar with 
cartatn additional variablaa partinant to tria — tha occurranca of 
limiting* tha axiatanca of tha vortax ring atata on ona or mora propal- 
lara or totora, ate* 

Following aach trim atata (whan ao raquaatad) is m liating of tha 
atability and raaponaa darlvativa array#* 

Whan a tima hiatory haa baan raquaatad, additional liatinga of vari- 
ablaa follow, ona (or aach print Intarval baginning with tha initial 
trim condition at t • 0* Tha variablaa ara tha aaaa aa thoaa printad 
out tor tha varioua trim condition# callad for* Thia continuaa until 
tha apacifiad ainulation tima occur#. 
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OVEKVIBH 

The HYBRDS programs ver« written with future expanalon possibilities 
In 'alnd, depending upon the nature of the application. The aajorlty of 
tl;c extensions or elaborations are quite straightforward* aeanxng that 
the subroutines Involved and the Interfaces with che rsaalnlng prograa 
are few in nuaber. This section briefly describes the nature of the 
changes required for the elaborations envisioned at the time of writing. 

lORCK AND MDMEHT GBAKKS 

The equations for the forces and aoments acting on the hull assembly 
and the payload can be altered cr elaborated relatively easily when 
there are no changes In the independent state variables which are inte- 
grated forward in time. The guideline here Is to consider the operation 
of the trimming routines. If the forces and moments are quite sensitive 
to changes In the state variables, the trim routine can run Into trou- 
ble, e.g., with the landing gear In ULAMOR and with the payload suspen- 
sion cables In HLAPAY. 

Multipoint Mooring 

The present simulation assumes an ideal translational constraint at 
a single point and can trim the vehicle with one or more landing gears 
touching the ground. Additional restraining cables, modeled as spring- 
damper combinations, can be added relatively easily If one Is careful to 
consider how the trim routine will operate. 
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PsaylMd Gtoond Contact 


When ground contact forces are added to the payload It will beccae 
possible to simulate the changing flight conditions associated with load 
pickup and placement* The pickup case Is more coapllcated as It re- 
quires the trimming sequence to be altered. The payload attitude must 
be trimmed on the ground with ground contact forces providing the con- 
trols for the trim routine's operation* 

iDDITlOHAL STATE VASTAnT.itc 

The simulation already has provision for additional state variables 
In a number of areas* Taking advantage of these provisions Is straight- 
forward to the extent that the trimming function Is left unchanged* 

Flight Control System 

The present flight control system has six state variables requiring 
integration which are associated with the Integral equalization (and 
control surface autotrlm) of the system* Additional states involving 
filtering, l*e*, sensor filters, actuation lag filters, etc*, can be 
added — spare elements are available in the vector of state variable 
rate of change for the HLASIM and RLAPAY programs* Changes In the loop 
structure apart from adding additional states are also possible* 

Propulsion ^stem 

The dynamics of a governor-controlled turbine-powered propulsion 
system can be added easily If the mathematical model trim operating 
point la completely specified by the input data* Consider a three-state 
propulsion system model for each LPU where rotor speed, primary turbine 
speed, and fuel flow are the three states* If simple enough, the rotor 
speed Implies the trim values of the remaining variables without the 
Iterative operation of a trimmer being required* The propulsion system 
dynamics can be thought of as additional flight control loops where the 
trim control deflections, rotor speeds, are specified at the outset* 
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LPD Orientation Belatlve to the Bill 


In the early stages of the HYBRDS development, one of the mathema- 
tical model features was controllable LPU angles, as In, for example, 
the XV-IS Tilt Rotor Aircraft* The structure of the nultlbcdy equations 
of motion Is such that three glmbal angle accelerations are potential 
control variables at each LPU attach point* In fact, the angular accel- 
erations and velocities are still Integrated In the program* A program 
statement sets the velocity Integrand to zero each time step to avoid 
drift in the glmbal angles* 

If the mathematical model for LPU glmbal angle control Is such that 
glmbal angle acceleration (rather than torque) is the controlled vari- 
able, mcdest changes to the program will restore these degrees of free- 
dom* A further requirement Is that the glmbal acceleration and velocity 
both be zero and the glmbal angle be specified at trim to avoid elabor- 
ating the trim routine operation* Such a modeling is quite adequate for 
Irreversible controls such as screw Jacks* The program will Integrate 
the accelerations and velocities and compute the constraining torque as 
a dependent variable* 

Trim Control with Ballonets 

The HLA hull and tall assembly is currently modeled as having con- 
stant Inertial properties, the latter Including the so-called apparent 
mass terms which arise In the nonsteady and acceleration-dependent terms 
in the aerodynamic forces* However, alrshlp-pecullar properties asso- 
ciated with the expansion, contraction, and distribution of buoyant gas 
within the envelope are not present* 

Adding these effects requires that the mass matrix be Inverted each 
time step* Currently, this Is done once as a part of the input data 
processing* In addition, center-of-gravlty offset and center-of- 
buoyancy offset (from the nominal in each case) terms need to be added 
to the eauations of motion* Finally, the trim routine logic requires 
modification of ballonet air distribution as well as rotor thrust •— 
both used for control of pitching moment by the flight control system. 
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HYBRDS has two post-processors • one for formatting plotting files 
and the other for computing additional gust response derivatx^re matri- 
ces • There are no pre-processor routines* Clearly anything in the 
extensive Input data set Is a potential candidate for pre-processing; 
and there are obvious possibilities for post-processing. An example of 
each Is cited below. 

Gust Inputs 

The simulation accepts the equivalent of prerecorded data for the 
disturbance Inputs. That the data be prerecorded Is a requirement of 
the error-correcting integration routine which alters step size as 
required to satisfy Its error criterion. There are several sets of 
these data corresponding to the four (or five in HLAPAY) gust "sources" 
of the model. Tb^ nature of the data (amplltudet random versus deter- 
ministic characteristics, correlation among the sources, etc.) is user 
specified* Given the appropriate models of these disturbances, one 
could create a pre-processing routine to generate the required data 
flies based on user-specified parameters. 

load Distribution 

The simulation generates as a part of Its output the Internal forces 
and moments acting at a number of points In the structure as well as a 
number of external force summations. Internal forces Include: 

• LPU attach point constraint forces and moments 

• Payload suspension cable tensions 

• Landing gear compressive loads 

• Mooring attach point loads 

• Rotor and propeller forces and torques 
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External force summations Include the forces and moments acting at the 
several areodynamlc centers In the system — the tall, the hull center 
of volume, the payload and LPU reference centers* 

The simulation does not perform a piecewise Integration of the 
forces acting along the length of the hull. However, within the accu- 
racy of the assumptions making up the aerodynamic force model, one can 
"work backwards" to Infer such data from the one thousand or so data 
elements available at each time step. From these data, hull bending 
moment and related loads can be Inferred. This Is a candidate for a 
post-processing routine. 

One could, of course, add this processing to the simulation Itself. 
With additional changes, this would allow computation of transfer func- 
tions relating, e.g., bending loads, to control Inputs and disturbances 
for linear systems analyses. 
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